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ABSTRACT 

Current  coronal  mass  ejection  (CME)  models  set  their  lower  boundary  to  be  in  the  lower  corona. 
They  do  not  calculate  accurately  the  transfer  of  free  magnetic  energy  from  the  convection  zone  to 
the  magnetically  dominated  corona  because  they  model  the  effects  of  flux  emergence  using  kine¬ 
matic  boundary  conditions  or  simply  assume  the  appearance  of  flux  at  these  heights.  We  test  the 
importance  of  including  dynamical  flux  emergence  in  CME  modeling  by  simulating,  in  2.5D,  the 
emergence  of  sub-surface  flux  tubes  into  different  coronal  magnetic  field  configurations.  We  inves¬ 
tigate  how  much  free  magnetic  energy,  in  the  form  of  shear  magnetic  field,  is  transported  from  the 
convection  zone  to  the  corona,  and  whether  dynamical  flux  emergence  can  drive  CMEs.  We  find 
that  multiple  coronal  flux  ropes  can  be  formed  during  flux  emergence,  and  although  they  carry  some 
shear  field  into  the  corona,  the  majority  of  shear  field  is  confined  to  the  lower  atmosphere.  Less 
than  10%  of  the  magnetic  energy  in  the  corona  is  in  the  shear  field,  and  this,  combined  with  the  fact 
that  the  coronal  flux  ropes  bring  up  significant  dense  material,  means  that  they  do  not  erupt.  Our 
results  have  significant  implications  for  all  CME  models  which  rely  on  the  transfer  of  free  magnetic 
energy  from  the  lower  atmosphere  into  the  corona  but  which  do  not  explicitly  model  this  transfer. 
Such  studies  of  flux  emergence  and  CMEs  are  timely,  as  we  have  new  capabilities  to  observe  this 
with  Hinode  and  SDO,  and  therefore  to  test  the  models  against  observations. 

Subject  headings:  CMEs,  Flux  Emergence,  MHD 


1.  INTRODUCTION 


1.1.  CME  modeling 


Coronal  mass  ejections  (CMEs)  and  eruptive  flares  are  the  most  energetic  manifestations  of  solar  ac- 
tivity,  with  a  typical  CME  accelerating  1016  g  of  coronal  plasma  to  speeds  sometimes  exceeding  1000  km/s 
( Gosling  et  al.  1976;  Hildner  197 3).  It  is  now  generally  accepted  that  these  giant  eruptions  of  solar  plasma  and 
field  are  due  to  the  explosive  release  of  magnetic  energy  stored  in  the  corona  prior  to  the  event.  Furthermore, 
the  free  energy  for  CMEs /flares  is  believed  to  be  stored  in  the  strongly  sheared  flux  of  filament  channels  (see, 
e.g,  reviews  by  Forbes  (2QQol):  Klimchuk  (20011:  Linton  fc  Moldwin  ( 2009^) ) .  This  flux  is  sheared  in  that  that  the 
field  lines  are  mainly  parallel  to  the  photospheric  neutral  line,  rather  than  perpendicular  as  would  be  expected 
for  a  potential  field.  All  CMEs/flares  are  associated  with  filament  channels  and  these  are  the  only  locations 
in  the  corona  where  the  magnetic  field  exhibits  strong  non-potentiality.  The  basic  picture  of  a  CME  is  that, 
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prior  to  the  event,  there  is  a  force  balance  in  the  corona  between  the  upward  magnetic  pressure  of  the  sheared 
filament  channel  held  and  the  downward  tension  of  overlying  unsheared  held.  As  the  shear  builds  up  or  the 
overlying  held  evolves,  this  force  balance  eventually  disrupts  producing  an  explosive  outward  expansion  of  the 
hlament  channel  flux  and  some  of  the  overlying  flux. 

It  is  evident  from  this  discussion  that  the  two  central  issues  in  understanding  CMEs/eruptive  hares  are 
the  process  by  which  magnetic  shear  forms  in  hlament  channels,  and  the  mechanism  that  disrupts  the  force 
balance.  These  two  issues  have  been  at  the  heart  of  CME  and  hare  research  for  decades  and  are  the  features 
that  distinguish  the  various  CME  models.  The  hrst  issue  of  the  shear  formation  process  is  especially  critical, 
because  the  disruption  mechanism  is  almost  certain  to  depend  on  the  topology  of  the  hlament  channel  held, 
which  must  be  determined  by  the  formation  process.  For  example,  some  models  invoke  kink-like  instabilities, 
which  requires  the  formation  of  magnetic  twist  as  well  as  shear  in  the  hlament  channel  held.  We  emphasize, 
however,  that  all  the  models  and  all  the  observations  agree  that  the  bulk  of  the  magnetic  free  energy  is  in  the 
form  of  shear.  Note  also  that  the  issue  of  shear  formation  has  broad  implication  for  a  number  of  important  solar 
physics  problems,  in  particular  the  origin  of  prominences  and  filaments. 


There  are  only  two  general  processes  that  can  produce  the  observed  shear  in  hlament  channels:  either  the 
held  emerges  from  below  in  a  non-potential,  sheared  state,  or  the  photospheric  motions  shear  the  held  after  it 
emerges  into  the  corona.  Of  course,  both  processes  must  be  present  to  some  extent  in  the  real  corona.  The 
latter  process  is  likely  to  dominate  in  non-active  region  hlament  channels,  such  as  those  associated  with  high- 
latitude  quiescent  hlaments,  because  these  form  well  after  the  flux  has  emerged.  On  the  other  hand,  in  active 
regions,  especially  in  the  strong  complex  regions  that  are  the  source  of  fast  CMEs,  the  hlament  channels  form 
with  the  active  region  and,  hence,  the  shear  must  emerge  with  the  flux.  Given  that  the  shear  formation  is  the 
fundamental  driver  of  the  eruption,  it  is  clear  that  flux  emergence  must  be  explicitly  included  in  any  model  for 
fast  CMEs  in  order  for  the  model  to  be  physically  rigorous. 


One  such  model,  which  generates  fast  CMEs  from  sheared  hlaments,  but  which  should  be  modihed  to 
include  the  creation  of  these  sheared  hlaments  via  flux  emergence,  is  the  magnetic  breakout  model  (Antipchos 
1998:  Antiochos  et  al.lfl999h.  This  model  exploits  the  finite  resistivity  in  the  solar  atmosphere  to  disrupt  the  pre- 


CME  force  balance  by  generating  reconnection  external  to  the  hlament  channel.  This  reconnection  reconhgures 
the  magnetic  held  and  allows  the  eruption  of  the  highly  sheared  structure.  The  breakout  model  thus  requires 
a  multipolar  flux  distribution  with  a  separatrix  where  reconnection  can  occur.  Numerical  simulations  have 
investigated  this  model  in  both  2.5D,  where  the  domain  is  axisymmetric  and  two  dimensional  but  where  all 
three  components  of  vector  variables  evolve  ( Lynch  et  alJl2004  [  MaeNeice_et_al|2004  [  DeVore  &;  Antiochosll2005 ) 


and  3D,  where  the  domain  is  three  dimensional  (Lynch  et  al.  2008:  DeVore  fe  AntiochoJ20081.  The  basic  model 


uses  a  current-free  quadrupole  magnetic  held  configuration,  consisting  of  a  central  active  region  dipole  embedded 
in  an  anti-parallel  global  dipole  held,  with  a  null  point  separating  the  two  systems.  This  is  shown  in  Figure 
ma),  taken  from  iLvnch  et  all  |2004).  Figure  (Ha)  shows,  for  this  experiment,  the  vertical  force  balance  in  the 
central  arcade  along  the  latitude  =  0°  symmetry  line  for  the  initial  equilibrium.  The  gravitational  (pg)  and 
gas  pressure  (—VP)  forces  balance,  as  do  the  magnetic  tension  (B.VB/47t)  and  magnetic  pressure  ( ~'VB2/8ir ). 
Shear  hows  are  applied  at  the  lower  boundary,  within  the  central  arcade,  to  generate  magnetic  held  parallel  to 
the  active  region  neutral  line  (in  our  2.5D  simulations  the  shear  field  will  be  the  magnetic  held  in  the  ignorable 
direction).  Figure [U(b)  shows  the  vertical  forces  after  a  short  period  of  shearing.  The  presence  of  shear  held 
exerts  an  upward  magnetic  pressure  on  the  overlying  held  and  stretches  the  original  null  point  into  a  current 
sheet,  as  shown  in  Figures  [lib)  and  (c).  Reconnection  occurs  when  the  current  sheet  collapses  below  grid  scale 
and  overlying  held  is  removed  by  this  reconnection,  which  allows  further  expansion  of  the  central  arcade.  The 
outward  expansion  of  the  sheared  arcade  occurs  exponentially  as  more  and  more  overlying  held  is  reconnected. 
Flare-like  reconnection  underneath  the  sheared  dipole  creates  a  disconnected  hux  rope  which  can  then  escape, 
as  seen  in  Figures [Hd)-(f),  at  speeds  comparable  to  observations  of  a  fast  CME. 


There  are  numerous  other  CME  models  which  either  assume  the  appearance  of  sheared  magnetic  hux  at  the 
low  corona,  or  use  driving  velocities  to  create  sheared  structures  from  coronal  equilibria.  Chen  fc  Shibatai  (200(il 
and  Dubev  et  al.  (20061)  simulated  a  2.5D  CME  model  where  a  pre-existing  coronal  hux  rope,  constrained  by 
overlying  held,  is  allowed  to  escape  when  kinematically  imposed  emerging  hux  causes  reconnection  in  a  hlament 
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channel  beneath  the  flux  rope.  They  also  found  that  CME  eruptions  occur  when  kinematic  flux  emergence 
is  used  to  drive  reconnection  outside  the  filament  channel.  The  flux  cancellation  model  ( Amari  et  al.  200d 


2003, 12007)  is  another  model  based  on  reconnection,  but  this  model  uses  both  boundary  imposed  shearing  and 
the  cancellation  of  magnetic  flux  at  the  neutral  line  of  filament  sites  to  initiate  a  CME.  Fan  fc  Gibson  (120071) 
developed  a  CME  model  where  the  lower  corona  is  driven  by  the  kinematic  emergence  of  a  sheared  flux  rope 
into  a  pre-existing  coronal  field.  Here,  loss  of  equilibrium  via  the  kink  or  torus  instability  leads  to  an  eruption, 
with  the  strength  of  the  overlying  field  being  an  important  constraint  on  the  eruption.  The  breakout  model  has 
also  been  extended  to  cover  different  magnetic  field  configurations,  a,nd  also  to  include  a  solar  wind,  in  both 


2.5D  and  3D  simulations  (van  der  Holst  et  al.ll2007t  IZuccarello  et  al.ll2008t  ISoenen  et  al.ll2009f).  Other  CME 


models  exist  which  do  not  directly  rely  on  flux  emergence,  such  as  those  based  on  the  loss  of  equilibrium  of  a 
coronal  flux  rope  (Titov  fc  Demoulin  19991:  Roussev  et  al.  2003  L  Klien^&Torok  2006),  and  simulations  of  the 
rise  of  a  pre-formed  coronal  flux  rope  into  the  solar  wind  by  Zhang  fc  Wu  ( 20091) .  Although  not  relying  on  flux 
emergence  directly,  these  models  do  assume  that  the  flux  ropes  have  already  been  formed  in  the  corona  and 
reached  an  equilibrium  with  the  background  coronal  field,  a  state  which  is  most  likely  caused  by  the  emergence 
of  magnetic  field  from  beneath  the  surface. 

All  these  models,  constrained  by  the  need  to  extend  the  simulation  domain  to  at  least  a  few  solar  radii, 
do  not  model  the  lower  solar  atmosphere.  The  lower  boundary  of  these  simulations  has  a  typical  density  of 
3  x  10-16  g/crn  and  a  typical  plasma- /3  <  0.2  where 


P  = 


P 


B2/8n' 


(1) 


Based  on  the  VAL  model  of  the  solar  atmosphere  (jVernazza  et  al.  119811)  the  photosphere,  in  contrast,  has  a 
density  of  around  3  x  10~7  g/cm3  and  /?  >  1.  These  CME  initiation  models  simply  assume  the  appearance  of 
flux  at  the  low  density,  low  fi  corona  and  do  not  self-consistently  calculate  a  process  for  the  flux  emergence. 
These  simulations  do  not,  therefore,  address  the  critical  question  of  whether  and  how  newly  emerging,  sheared 
magnetic  flux  can  rise  from  its  origins  in  the  high  /?  convection  zone  to  the  low  corona  where  it  is  required  to 
drive  CME  models.  This  is  an  important  question  to  address  now,  as  Hinode  and  SDO  scientists  are  currently 
making  detailed  observations  of  flux  emergence  with  vector  magnetograms,  which  will  allow  the  comparison  of 
observations  and  state  of  the  art  CME  models. 


1.2.  Flux  Emergence 


The  current  state  of  the  theory  of  the  formation  of  active  regions  suggests  that  dynamo  actions  in  the  solar 
interior  create  twisted,  buoyant  flux  tubes  which  rise  through  the  convection  zone,  intersect  the  photosphere 


and  create  the  observed  sunspots  and  bipolar  active  regions  (IParker 


1955 


Zirin  197(ll:  Schuessler  1979:  Parker 


1988:  Zwaan  1987:  Emonet  fc  Moreno-Insertis  19981:  Fisher  et  al.1  2000).  The  subsequent  expansion  of  these 
flux  tubes  into  the  corona  has  been  a  subject  of  vigorous  research  for  the  past  30  years.  The  most  likely 
driver  of  flux  emergence  into  the  lower  corona  is  an  ideal  MHD  instability  driven  by  magnetic  buoyancy  (the 


magnetic  buoyancy  instability),  a  particular  mode  of  which  was  first  suggested  by  IShibata  et  al 


2D  simulations  using  a  crude  model  atmosphere  and  magnetic  field  configuration.  They  found  that  the  rise  of 
magnetic  flux  into  the  corona  due  to  this  instability  matched  observed  rise  velocities  of  magnetic  flux. 

The  process  of  flux  emergence  has  been  studied  more  recently  in  3D  simulations  with  a  more  detailed 
model  atmosphere  that  includes  a  convection  zone,  isothermal  photosphere/chromosphere,  transition  region,  and 
isothermal  corona,  with  a  buoyant  twisted  flux  tube  embedded  in  the  convection  zone  ( Fanll2001 ;  Manchester  et  al 
2004:  Archontis  et  al.  12004  iFan  12009).  The  reader  is  directed  to  a  comprehensive  review  of  this  subject  by 
Archontis  (120081) .  These  simulations  show  that  as  the  flux  tube  intersects  the  photosphere,  bipolar  active  re¬ 
gions  are  formed  which  exhibit  shear  flows  along  their  neutral  line.  As  the  tube  continues  to  emerge  into  the 
corona,  signroid-like  structures  are  formed  consistent  with  observations,  with  magnetic  field  being  created  par¬ 
allel  to  the  neutral  line  (shear  field).  This  suggests  that  flux  emergence  may  be  able  to  provide  magnetic  shear 
field  sufficient  to  drive  CMEs,  but  it  is  unclear  how  this  magnetic  shear  emerges  into  the  corona,  as  the  axis  of 
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the  emerging  flux  tubes  remains  rooted  near  the  surface. 


Manchester  et  ah  02004)  showed  that  shearing  motions  associated  with  the  emergence  are  able  to  create 
a  sheared  coronal  flux  rope  which  rises  out  of  the  simulation  domain,  although  the  authors  do  not  associate 
this  with  an  eruption  because  the  rise  speed  is  limited  by  the  amount  of  available  ’axial  flux’  to  drive  an 
eruption.  Note  that  in  this  paper  we  use  the  term  flux  tube  to  refer  to  the  sub-surface  tubes  which  we  initiate 
our  simulations  with,  and  flux  rope  to  refer  to  the  twisted  structures  formed  in  the  corona.  3D  simulations 
(Fan  2009)  have  shown  that  without  a  pre-existing  coronal  field,  the  rise  of  coronal  flux  ropes  created  during 
flux  emergence  is  constrained  by  the  expanding  emerging  field  of  the  original  flux  tube,  or  envelope  field.  In  a 
similar  manner,  the  work  oflArchontis  fe  Hoodl  (1200811  andlMacTaggart  &;  Hoodl  (|2009af)  showed  that  flux  ropes 
formed  from  multiple  flux  tube  emergence  are  ultimately  trapped  by  the  envelope  fields.  For  an  eruption  to 
occur,  the  overlying  tension  of  the  envelope  field  must  be  removed. 


Archontis  fe  Torok  ([20081)  added  a  pre-existing  horizontal  coronal  field  to  their  flux  emergence  simulations. 
They  found  that  a  coronal  flux  rope  could  be  formed  by  inflows  and  horizontal  shearing  motions  which  drive 
reconnection  within  the  expanding  volume  of  the  flux  tube.  They  also  found  that  reconnection  between  the 
envelope  field  and  the  pre-existing  field  reduced  the  overlying  tension  of  the  envelope  field  and  allowed  the 
flux  rope  to  erupt.  MacTaggart  fc  Hoodl  (2009b)  found  a  similar  result,  but  used  a  toroidal  sub-surface  flux 
tube  which  was  already  arched  towards  the  surface,  which  allowed  the  flux  tube  axis  to  emerge  bodily  into 
the  atmosphere.  Archontis  &  Hood  (20101)  showed  that  the  likelihood  of  a  coronal  flux  rope  escaping  when 
reconnection  occurs  between  the  ambient  coronal  field  and  the  envelope  field  depends  strongly  on  the  strength 
of  the  ambient  field. 


Flux  emergence  simulations  such  as  these  occur  over  timescales  of  1000  s,  much  quicker  than  timescales 
of  CME  models,  and  do  not  extend  high  enough  in  the  corona  to  allow  a  full  comparison  with  CME  models 
and  observations.  The  development  of  coupled  models  which  span  the  timescales  and  spatial  scales  of  both  flux 
emergence  and  CME  eruption  is  very  important  to  fully  understand  the  initiation  of  these  explosive  events. 

In  this  work  we  investigate  whether  magnetic  flux  emergence,  in  a  2.5D  cartesian  geometry,  can  provide 
enough  magnetic  shear  energy  to  drive  CME  initiation,  with  focus  on  the  initiation  mechanism  of  the  magnetic 
breakout  model.  Figure  Oc)  shows  the  energy  in  the  shear  magnetic  field  (normalized  to  the  energy  in  the  total 
field)  in  a  typical  breakout  model  simulation  as  shearing  is  applied  at  the  surface.  Also  shown  is  the  total  kinetic 
energy  in  the  system.  As  the  inner  arcade  is  sheared  at  the  ft  =  0.2  lower  boundary  the  shear  energy  increases. 
When  the  magnetic  shear  energy  is  about  20%  of  the  total  magnetic  energy,  the  upward  magnetic  pressure  due  to 
the  shear  field  can  force  reconnection  at  the  null  point  and  hence  further  expansion,  resulting  in  a  rapid  increase 
in  kinetic  energy.  Hence  the  breakout  model  typically  needs  20%  of  the  magnetic  energy  above  the  f  =  0.2 
corona  to  be  in  the  shear  field.  Although  we  know  that  flux  emergence  can  create  the  shear  flows  and  sheared 
magnetic  fields  near  the  photosphere  which  are  similar  to  those  required  in  the  low  corona  by  CME  models, 
it  has  not  yet  been  calculated  how  much  shear  energy  the  emergence  of  magnetic  flux  can  actually  supply  to 
the  low  /3  corona.  This  issue  is  key  to  understanding  how  to  couple  dynamical  flux  emergence  simulations  and 
current  CME  simulations.  We  will  investigate  this  using  simulations  of  the  evolution  of  sub-surface  flux  tubes 
with  a  variety  of  tube  profiles  as  they  emerge  into  both  field-free  corona  and  corona  with  pre-existing  fields.  We 
will  also  investigate  whether  flux  emergence  can  drive  CME  initiation  in  these  configurations. 

This  paper  is  structured  as  follows:  §2  explains  the  model,  including  numerical  method,  atmospheric  and 
magnetic  models,  and  initialization.  §3  presents  results  from  2.5D  simulations.  Firstly  we  look  at  emergence 
into  a  field-free  corona,  then  we  look  at  emergence  into  a  quadrupole,  and  finally  emergence  into  a  dipole.  In 
all  cases  we  attempt  to  cover  a  realistic  range  of  field  profiles  and  coronal  field  strengths  to  test  the  robustness 
of  our  results.  §4  summarizes  our  findings  and  suggests  consequences  for  modeling  the  initiation  of  CMEs  by 
flux  emergence. 
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2.  MODEL  DESCRIPTION 


2.1.  Numerical  Method 


We  model  our  system  of  a  magnetized  plasma  with  a  polytropic  equation  of  state  (polytropic  index  7  =  5/3) 
using  the  resistive  MHD  equations.  These  are  solved  numerically  using  the  Lagrangian  remap  code  LareXd 
(lArber  et  alJl2001).  The  equations  solved  are  presented  below  in  Lagrangian  form,  using  Gaussian  units: 


=  — pV.v 

(2) 

1  „  1  j  A  B  1  „ 

(3) 

VP+  J  +g+  v.s 

p  p  c  p 

=  (B.V)v  —  B(V.v)  —  cV  A  {pj) 

(4) 

Pv7  ,  -2  ,  0  e-e0(p) 

—  V.V  TJJ  +  Qj  Sij 

(5) 

Dp 
Dt 
Dv 
~Dt 
DB 
~Dt 
De 
Dt. 

The  gas  density,  pressure,  and  internal  specific  energy  density  are  denoted  by  p,  P,  and  e  respectively  and 
are  defined  at  the  center  of  each  numerical  cell.  The  magnetic  field,  denoted  by  B,  is  defined  at  cell  faces, 
j  =  cV  A  B/4-7T  is  the  current  density,  and  c  is  the  speed  of  light.  The  velocity,  v,  is  defined  at  cell  vertices. 
This  staggered  grid  preserves  V.B  during  the  simulation.  The  gravitational  acceleration  is  denoted  by  g,  v  is 
the  viscosity,  set  to  2.6  x  104  g.cm_1s_1,  p  is  the  resistivity,  and  S  is  the  stress  tenor  which  has  components 
Sij  =  ij  —  g^jV.v),  with  Qj  =  +  fy2-)-  The  four  equation  are  closed  with  a  simple  equation  of  state 

representing  an  ideal  gas:  P  =  pRT ,  where  R  is  the  gas  constant. 

Equation [5j  minus  the  last  term,  (e  —  eo (p))/r,  is  similar  to  that  used  in  the  majority  of  flux  emergence 
simulations  to  date  ( Archontisl[2008 ) .  When  flux  tubes  emerge  and  expand  into  the  corona,  the  associated  cooling 
due  to  the  pressure  terms  can  lead  to  unrealistic  temperatures  in  the  corona,  as  shown  by  Leake  fe  Arbeil  (2006  b 
Leake  fc  Arben  (|200a)  also  showed  that  by  including  a  relaxation  term,  such  as  the  last  term  in  Equation  0 
this  unrealistic  cooling  could  be  avoided.  This  approach  is  a  simple  way  to  model  the  effects  of  terms  that  we 
are  currently  unable  to  formulate  numerical  equations  for  (such  as  chromospheric/coronal  heating),  and  terms 
such  as  thermal  conduction  and  radiative  losses  which  are  beyond  the  scope  of  this  particular  numerical  code. 
The  term  (e  —  eo(p))/r  relaxes  the  specific  internal  energy  density  back  to  its  initial  equilibrium  values  eo(p)  on 
a  timescale  r  which  scales  with  density: 


r(y)  =  23 


p{y ) 

p{y  =  0) 


-1.43 


s. 


(6) 


The  equations  are  solved  in  2.5D:  the  simulation  box  is  2D,  with  x  and  y  being  independent  variables  and 
z  being  ignorable,  but  all  three  components  of  the  vector  variables  are  evolved.  The  simulation  box  extends 
from  -3  Mm  to  90  Mm  in  the  vertical  direction  (y)  and  from  -45  Mm  to  45  Mm  in  the  horizontal  direction  (x). 
This  is  much  larger  than  most  flux  emergence  simulations  typically  cover  (Manchester  et  al.  2004:  Fan  2009h. 
but  is  designed  to  allow  expansion  of  the  flux  tubes  into  coronal  field  as  far  as  possible.  The  numerical  grid  is 
stretched  in  both  x  and  y  to  provide  better  resolution  in  areas  of  interest.  The  grid  spacing,  8,  is  smallest  at  the 
surface,  where  8  =  0.0375  Mm,  0.25  times  the  scale  height  at  the  surface,  and  largest  at  the  upper  boundary, 
where  8  =  0.15  Mm. 

In  the  majority  of  previous  flux  emergence  simulations  either  ideal  MHD  is  used  {p  =  0),  or  p  is  chosen 
to  be  a  constant  such  that  the  diffusion  due  to  this  explicit  resistivity  is  greater  than  the  numerical  diffusion 
in  the  code.  Theoretical  evidence  suggests  that  when  the  electron  fluid  flow  speed,  ve,  exceeds  the  phase  speed 
of  the  ion-acoustic  model,  cla .  then  ion- acoustic  turbulence  has  a  strong  effect  on  current  sheet  development 


([Bvchenkov  et  ajJ  1 


In  this  case  an  effective  formula  for  anomalous  resistivity  would  be 


p  =  po  max  (  0, - 1 

C-ia 


(7) 
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which  we  can  rewrite,  using  ve  ~  |j|/rae,  as 


77  =  770  max  0 


-  1 


Jcrit 


(8) 


where  jcrn  =  necia.  We  choose  770  =  1.6  x  10-11  s,  and  jcru  =  2.07  x  107  statampere.cm-2,  and  this  anomalous 
diffusion  exceeds  the  numerical  diffusion  in  the  code.  This  approach  has  been  used  in  MHD  simulations  of  kink 
instabilities  of  coronal  loops  (Arber  et  al.  1999k  Gerrard  et  al.  2001). 


2.2.  Initial  Conditions 


2.2.1.  Background  Stratification 


The  initial  stratification  is  a  simple  ID  model  of  the  temperature  profile  of  the  Sun  which  includes  the  upper 
3  Mm  of  the  convection  zone,  a  photosphere/chromosphere,  transition  zone,  and  corona.  The  temperature  profile 
is  given  by 


T(y) 


\  Tph 

Isl 

R(m+ 1)  »> 

[  Tph  + 

(Tcor— Tph) 

2 

tanh  ( y  ytr)  +  1 

\  ^tr  ) 

y  <  0 
y>  0 


(9) 


The  convection  zone  profile  (y  <  0)  is  a  linear  polytrope  which  is  marginally  unstable  to  convection,  with 
771  =  t-A|-  being  the  adiabatic  index  for  a  polytrope.  The  temperature  in  the  photosphere  and  chromosphere  is 
assumed  to  be  constant,  Tph  =  5700  K,  as  is  the  temperature  in  the  corona,  Tcor  =  8.6  x  105  K.  The  height 
of  the  transition  region  is  ytr  =  3.75  Mm,  and  its  width  is  Wtr  =  0.75  Mm.  The  density  and  gas  pressure  are 
specified  by  initially  assuming  hydrostatic  equilibrium.  The  resulting  stratification  is  shown  in  Figure  (SJa). 


2.2.2.  Magnetic  Flux  Tube 


We  insert  a  cylindrical  magnetic  flux  tube  into  the  model  convection  zone  at  a  height  of  yt  =  —1.8  Mm. 
The  axial  field  for  all  the  tubes  is  given  by 

Bz  =  B0e~r2/a2  (10) 


where  r  =  y/x2  +  {y  —  yt )2  is  the  radial  distance  from  the  center  of  the  tube,  B0  =  7800  G  is  the  axial  field 
strength  at  the  center  (r  =  0),  and  a  =  0.3  Mm  is  the  width  of  the  tube.  This  axial  field  will  equate  to  the  shear 
field  required  by  the  breakout  model.  We  then  choose  three  different  tubes  by  choosing  three  different  profiles 
for  the  twist  magnetic  field  ( Bg  in  the  cylindrical  coordinate  system,  Bx  and  By  in  cartesian).  A  minimum 
amount  of  twist  is  required  to  play  the  role  of  surface  tension  and  keep  the  flux  tube  coherent  as  it  rises  in  the 
convection  zone  (Emonet  fe  Moreno-Insertisi  1998).  We  denote  the  twist  by  0,  defined  as 


0  = 


Be 

rBz 


The  three  twist  profiles  we  use  are  given  by 


Tube  1  :  0i(r)  =  c 

Tube  2  :  02  (r)  = 

Tube  3  :  03(r)  =  ce 

where  c  =  1/a,  Ci  =  0.15  Mm  and  c2  =  1.5  Mm. 


(1 -*-■'■*) 

-r2/4 


(11) 

(12) 

(13) 

(14) 


Tube  1,  the  constant  twist  tube,  is  similar  to  tubes  which  have  been  studied  extensively  in  both  2D  and 


3D  (Fan  2001:  Manchester  et  al.  2004:  Arber  et  al.  2007:  Leake  fe  Arber  2006:  Fan  20091).  This  initial  condition 
has  been  shown  to  create  a  buoyant  flux  rope  that  can  not  only  emerge  into  the  solar  atmosphere,  but  interact 
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with  a  simple,  horizontal  pre-existing  coronal  field  and  create  isolated  flux  ropes  in  the  corona,  and  is  thus  an 
ideal  choice  for  our  study. 


Tube  2  (increasing  twist)  and  Tube  3  (decreasing  twist),  are  similar  to  two  tubes  investigated  bylMurray  &  Hood; 


(120081)  and  have  twist  profiles  which  increase  and  decrease  with  radius,  respectively.  In  their  simulations  inves¬ 
tigating  the  effect  of  non-constant  twist,  the  tops  of  these  two  flux  tubes  emerged  higher  into  the  atmosphere 
than  tubes  with  other  twist  profiles  and  are  therefore  the  most  relevant  to  this  study. 


As  in  Murray  fe  Hood  ( 2008h .  we  can  group  the  three  tubes  into  low  magnetic  tension  and  high  magnetic 


tension  cases.  The  magnetic  tension  for  these  tubes  is  given  by 


(B.V)B  _  (■ -Be2/r )  er 
47T  47T 


The  magnetic  field,  twist  (0),  and  tension  profiles  for  the  three  tubes  are  shown  in  Figure |4]  Tubes  1  and  3  have 
similar  twist  and  tension  profiles.  Tube  2  has  lower  tension  near  its  axis  due  to  the  lower  magnitude  of  Bg  close 
to  the  center  of  the  tube.  From  their  3D  simulations,  Murray  fe  Hood  (2008lf  concluded  that  the  twist  profile 
was  not  as  important  a  factor  in  the  emergence  of  the  flux  tubes  as  the  axial  field  strength,  which  affected  the 
rate  of  emergence  but  not  the  extent  of  the  emergence.  The  purpose  of  choosing  these  three  tubes  in  our  study 
is  to  cover  a  range  of  flux  tube  behavior  whilst  keeping  the  potential  amount  of  shear  energy  provided  by  flux 
emergence  optimal. 


The  plasma  f3  at  the  center  of  the  flux  tubes  is  approximately  10,  and  the  flux  tubes  are  all  non-force- 
free.  To  initiate  a  buoyant  flux  tube  we  perturb  the  background  gas  pressure  by  an  amount  p\  such  that 
(Vpi)r  =  (j  A  B)r,  which  can  be  integrated  to  give 


Pi(r) 


AW  +  f_2L* 

Sn  Jr  (47t  r) 


(16) 


so  that  the  tube  is  in  radial  force  balance.  Assuming  that  the  flux  tube  is  in  thermal  equilibrium  with  its 
surroundings  makes  the  tube  less  dense  that  the  surrounding  plasma  and  initiates  its  buoyant  rise  to  the 
surface. 


We  will  investigate  the  emergence  of  these  tubes  first  into  a  field-free  corona  and  then  into  different 
quadrupole  and  dipole  coronal  field  configurations  to  test  how  efficient  flux  emergence  is  at  providing  mag¬ 
netic  shear  energy  to  the  corona. 


2.2.3.  Quadrupole  Background  Magnetic  Field 


After  investigating  the  emergence  of  flux  tubes  into  a  field  free  corona,  we  investigate  the  interaction 
of  coronal  field  and  emerging  flux  within  the  context  of  magnetic  breakout.  To  do  this  we  first  impose  a 
quadrupole  coronal  field  above  the  magnetic  flux  tube.  We  construct  the  quadrupole  as  in  Karpen  et  al.  ( 1996lf 
by  first  defining  the  magnetic  vector  potential  along  a  source  surface  at  y  =  —7.5 Mm  (which  is  outside  the 
computational  domain) : 


Az{ x,  -7.5  Mm)  =  qd  ■ 


f  x2 

1  -  ( X/Xaf 

l  0- 

0  <  |x|  <  xa 
Id  >  xa 


(17) 


with  xa  =  21  Mm  being  the  horizontal  extent  of  vertical  qaudrupole  field  at  this  height.  The  vector  potential 
in  the  interior  of  the  domain  is  then  determined  by  integration  of  Laplace’s  equation  with  Equation  [l7]  as  the 
lower  boundary  condition  and  Az  =  0  at  y  =  oo.  The  quadrupole  consists  of  an  inner  dipole  and  an  overlying 
dipole  of  opposite  orientation,  separated  by  a  null  point  at  y  ~  3  Mm.  The  height  of  the  null  is  chosen  to  allow 
for  the  expansion  of  the  flux  tube  as  it  emerges  into  the  lower  corona,  during  which  the  cross-section  of  the  tube 
increases  by  a  factor  of  10  (jArber  et  al .112007).  Thus  the  0.3  Mm  flux  tube  will  stay  within  the  central  arcade 
of  the  quadrupole  when  it  emerges.  The  magnetic  configuration  of  the  sub-surface  flux  tube  and  overlying 


quadrupole  field  is  shown  in  Figure  [5]  For  the  breakout  mechanism,  most  of  the  horizontal  quadrupole  flux 
above  the  null  is  removed  by  reconnection  with  horizontal  flux  below  the  null.  We  therefore  vary  the  strength 
of  the  quadrupole  so  that  the  horizontal  flux  above  the  null  point,  along  the  x  =  0  line,  is  a  certain  fraction 
of  the  horizontal  flux  contained  between  the  center  of  the  sub-surface  flux  tube  and  the  null  point.  We  choose 
qd  to  be  [11.7,  5.85,  3.9]  xlO7  Gem,  which  give  quadrupole  surface  strengths  at  x  =  xa/2  of  [197,  98,  66]  G 
respectively.  The  horizontal  flux  above  the  null  point  is  then  [1,  0.5,  0.3]  times  the  horizontal  flux  between  the 
flux  tube  and  the  null  point.  This  allows  us  to  control  the  amount  of  reconnection  at  the  null  point,  whilst 
still  being  able  to  fix  the  /3  at  certain  heights  above  the  surface.  The  curves  in  Figure  E|b)  shows  the  /3  profiles 
for  the  configurations  of  a  sub-surface  flux  tube  plus  three  different  strength  quadrupoles.  These  profiles  are 
generally  consistent  with  /3  profile  models  of  active  regions  developed  by  Gary  fe  Alexander  ( 1999^1  and  Gary 
Q2001)  which  take  into  account  Soft  X-ray  Telescope  limb  observations. 


2.2-4-  Dipole  BackGround  Magnetic  Field 


While  emerging  a  flux  tube  into  a  quadrupole  is  the  first  step  to  coupling  flux  emergence  and  the  breakout 
model,  we  also  investigate  a  more  self-consistent  configuration,  with  a  flux  tube  emerging  into  a  background 
dipole  field.  By  using  a  dipole  with  opposite  orientation  to  the  upper  half  of  the  flux  tube,  reconnection  between 
these  two  systems  creates  a  quadrupole  structure  self-consistently. 

The  dipole  field  is  represented  by  the  vector  potential  A  =  Azez  where 


Az  (x,  y)  =  d 


y  -ya 


(18) 


with  r\  =  ^x2  +  (y  —  yd)2  being  the  distance  from  the  source.  We  chose  yd  to  be  -15  Mm  so  that  the  flux 
tube  is  far  from  the  source  of  the  dipole.  To  cover  various  dipole  strengths  we  pick  a  range  of  d  =  [18,  9, 
4.5]  x  1028  Gem3  which  gives  a  magnetic  field  strength  at  x  =  0 ,y  =  0  of  [105,  53,  26]  G  respectively.  In  a  similar 
fashion  to  the  quadrupole  field,  we  choose  these  values  so  that  the  horizontal  flux  contained  between  the  axis 
of  the  flux  tube  and  the  separatrix  which  separates  the  flux  tube  and  the  dipole  is  a  certain  fraction  of  the 
horizontal  flux  above  this  separatrix.  With  decreasing  dipole  field  strength,  these  factors  are  [1,  0.5,  0.25].  The 
/3  profile  for  the  three  dipole  choices  overlying  a  sub-surface  flux  tube  are  shown  in  Figure  [3jb) .  These  choices 
of  dipole  strength  allow  for  a  range  in  the  0  profiles,  but  are  still  consistent  with  the  models  of  0  in  the  solar 
atmosphere  developed  by  Gary  &  Alexander!  (jl999)  and  Gary  (120011. 


3.  RESULTS 

3.1.  Initial  Evolution  in  the  Convection  Zone  and  Lower  Atmosphere:  Effect  of  Twist  Profiles 


The  evolution  of  the  flux  tubes  in  the  convection  zone  is  similar  to  previous  simulations  of  rising  flux  tubes 
(Emonet  fc  Moreno- Inserti3l998:  Magaral2001 ;  Leake  fc  Arbeillioooll .  The  initial  buoyancy  causes  the  flux  tube 
to  rise  to  the  surface,  during  which  time  the  cross-section  increases  and  flux  conservation  decreases  the  axial,  or 
shear,  field  strength.  As  the  flux  tube  meets  the  convectively  stable  photosphere  its  rise  is  halted  and  horizontal 
expansion  spreads  the  flux  tube  out  to  form  a  contact  layer  with  the  plasma  above.  The  density  above  this 
layer  is  higher  than  the  density  below  due  to  the  concentration  of  magnetic  field  and  the  layer  is  unstable  to  a 
Rayleigh-Taylor  like  instability  known  as  the  magnetic  buoyancy  instability  (Newcombl  1961:  Athav  &  Thomasl 
1961:  Acheson|l979).  The  stability  of  the  contact  layer  involves  a  competition  between  the  destabilizing  gradient 
in  the  magnetic  field  and  the  stabilizing  sub-adiabatic  temperature  gradient  ( Gilman!  1 9  7ol).  The  stabilizing  term 
is  dependent  on  the  local  0,  and  as  the  center  of  the  tube  rises  up  to  this  contact  layer,  the  local  3  falls  and 


the  instability  allows  the  upper  portion  of  the  flux  tube  to  expand  into  the  atmosphere,  as  in  Archontis  et  ai 

fcoQlh. 


As  can  be  seen  in  Figure  [He),  the  choice  of  twist  profile  affects  the  amount  of  tension  in  the  tubes.  The 
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increasing  twist  tube  (Tube  2)  has  relatively  less  tension  near  its  axis  than  the  constant  twist  tube  (Tube  1) 
and  the  decreasing  twist  tube  (Tube  3),  and  also  a  higher  buoyancy  at  its  center.  This  difference  distinguishes 
the  evolution  of  the  increasing  twist  tube  from  the  other  two,  and  from  now  on  we  refer  to  the  low  tension 
case  to  mean  Tube  2  and  the  high  tension  case  to  mean  Tubes  1  and  3.  In  Figure  [G]  we  can  see  how  the  choice 
of  twist  profile  affects  the  initial  expansion  into  the  field-free  corona.  The  tubes  act  nearly  identically  as  they 
rise  buoyantly  in  the  convection  zone  and  also  in  the  lower  atmosphere,  as  the  magnetic  buoyancy  instability 
develops  and  the  outer  fieldlines  extend  into  the  corona.  However,  the  two  cases  then  differ  in  how  the  center 
of  the  tube  reacts.  In  the  high  tension  cases,  the  lower  fieldlines  near  the  center  of  the  tube  remain  near  the 
surface,  shown  in  Figure  (Hid),  and  the  shear  field  is  concentrated  at  the  center,  creating  a  single  neutral  line.  In 
the  low  tension  case,  the  decreased  tension  allows  the  center  of  the  flux  tube  to  rise,  forming  a  crescent  shape, 
shown  in  Figure  [Gjh).  During  this  evolution,  the  mass  in  the  center  of  the  tube  drains  to  the  lowest  possible 
location  which  is  at  either  end  of  this  crescent  shape.  The  coupling  of  shear  field  and  density  in  this  2.5D 
simulation  means  that  the  shear  field  also  concentrates  in  these  regions.  Note  that  this  emerging  structure  now 
has  three  neutral  lines  at  the  photosphere. 


3.2.  Emergence  into  a  Field-Free  Corona:  Effect  of  Twist  Profiles 


The  center  of  the  high  tension  flux  tube  remains  at  the  surface  for  the  duration  of  the  simulation.  We 
follow  the  continued  evolution  of  the  low  tension  flux  tube  in  Figure  [3  which  shows  the  shear  field  and  the 
log  of  density,  along  with  magnetic  fieldlines.  The  deformation  of  the  center  of  the  tube  into  a  crescent  shape 
forms  a  current  sheet  as  field  of  opposite  sign  is  brought  together  on  either  side  of  this  crescent.  This  enhanced 
current  triggers  the  anomalous  resistivity  and  reconnection  sets  in.  As  can  be  seen  from  Figures  0  b)  and  (e), 
the  reconnection  sites  are  at  either  side  of  the  apex  of  the  current  sheet  arch.  Plasmoids  are  formed  which 
are  expelled  along  the  current  sheet.  The  half  width  of  the  current  sheet,  l,  is  approximately  0.15  Mm,  and 
a  typical  distance  between  the  plasmoids  is  5  Mm.  Thus  the  ratio  of  half  current  sheet  width  to  the  length 
between  islands,  A,  is  approximately  1/ A  =  0.03.  Based  on  Priest  &  Forbes]  (12000IK  the  wavelength  range  for  the 
tearing  mode  is 

1  /  DA1/4  l 

Sfe)  <  ^  <  1/(2-)  ~  0.16  (19) 

where  D  is  the  magnetic  diffusivity  (D  =  rjc2 /Air)  and  va  is  the  Alfven  speed.  Using  values  from  our  simulations 
gives  the  lower  limit  as  0.01.  Hence  the  current  sheet  dimensions,  1/ A  =  0.03,  are  compatible  with  a  tearing 
mode  instability.  Later  in  time,  plasmoids  are  formed  closer  to  the  apex  of  the  current  sheet  and  coalesce, 
creating  an  isolated  flux  rope  at  the  center.  Both  shear  field  and  mass  are  trapped  in  this  flux  rope,  but  the 
majority  of  the  shear  field  concentrates  at  either  end  of  the  current  sheet,  in  the  low-lying  dips  in  the  magnetic 
field. 


To  determine  what  happens  to  this  coronal  flux  rope  which  is  formed  from  the  original  center  of  the  low 
tension  flux  tube,  we  can  examine  the  forces  in  the  flux  rope.  Figure  [8]  shows  the  four  vertical  forces  associated 
with  gravity,  magnetic  tension,  magnetic  pressure,  and  gas  pressure  along  the  central,  x  =  0,  line  at  two 
different  times.  The  forces  are  shown  for  the  initial  sub-surface  flux  tubes  in  Figures  |5{a)  and  (b),  where  the 
initial  hydrostatic  equilibrium  can  be  seen  in  the  large  scale  trends  for  pg  and  —VP,  and  the  magnetic  forces 
can  be  seen  over  the  range  -2.5  Mm  to  -1.5  Mm.  The  outward  directed  magnetic  pressure  force,  due  primarily 
to  the  shear  field,  dominates  over  the  inward  directed  tension  force  of  the  twist  field.  The  gas  pressure  force 
is  directed  inwards  to  balance  this  excess  magnetic  pressure  force.  Figures  [5f c)  and  (d)  show  the  forces  in  the 
corona  after  the  flux  tube  has  emerged.  For  the  low  tension  case,  the  forces  in  the  newly  formed  coronal  flux 
rope,  the  center  of  which  is  at  y  =  7  Mm,  can  be  seen  in  Figure  Od).  The  presence  of  mass  can  be  detected  in 
the  small  dip  in  the  gravitational  force.  Here,  in  contrast  to  the  configuration  of  the  initial  flux  tube,  the  inward 
directed  tension  of  the  twist  field  dominates  over  the  outward  directed  magnetic  pressure  force.  Thus  the  gas 
pressure  force  is  directed  outward  here,  to  balance  the  excess  tension  force.  This  shows  that,  in  contrast  to  the 
breakout  model,  shear  field  does  not  create  a  strong  unbalanced  upward/outward  pressure  force,  and  so  there  is 
no  impetus  which  can  drive  a  CME  eruption.  This  shows  the  dramatic  effect  the  dense  lower  solar  atmosphere 
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can  have  on  the  role  of  shear  in  driving  CMEs. 


As  far  as  we  are  aware,  the  mechanism  shown  here  for  the  creation  of  a  coronal  flux  rope  from  an  emerging 
flux  tube  whose  cross  section  distorts  into  a  crescent  shape  has  not  been  seen  before.  We  expect  that  such 
coronal  flux  ropes  will  form  from  3D  flux  tubes  which  emerge  with  a  large  radius  of  curvature  along  their  axis. 
A  large  radius  of  curvature  implies  that  the  draining  of  mass  along  the  axis  will  not  exceed  the  draining  of 
mass  perpendicular  to  the  axis  seen  in  these  2.5D  simulations,  and  so  the  behavior  will  be  similar.  However, 
3D  simulations  of  similar  flux  tubes  by  Murray  &  Hood!  (2008)  do  not  show  this  behavior.  We  find  that  we 
still  see  the  formation  of  our  coronal  flux  ropes  when  we  use  the  same  parameters  (yt=- 0.17  Mm,  Ho=3900  G, 
a=0.425  Mm,  c  =  3  Mm-1  and  Ci=0.17  Mm)  as  one  of  the  low  tension  tubes  used  in  Murray  &  Hood  (120081) . 
3D  simulations  are  required  to  investigate  this  issue  further. 


Figure  [Ha)  shows  the  height  of  the  centers  of  the  three  flux  tubes  as  they  emerge  into  the  field- free  corona. 
The  center  (which  is  initially  the  axis)  is  defined  at  the  point  where  Bx \x-o  changes  sign.  The  centers  of  the 
high  tension  tubes  stay  at  the  photosphere,  while  the  center  of  the  low  tension  tube  rises  to  a  height  of  about 
8Mm  as  the  tube  is  deformed  into  a  crescent  shaped  current  sheet  and  the  resulting  coronal  flux  rope  rises.  The 
coronal  flux  rope  does  not  escape,  most  likely  due  to  the  presence  of  mass  accumulation  in  the  rope  and  due  to 
the  lack  of  magnetic  shear. 


Figures  [U(b)  and  (c)  show  the  normalized  magnetic  shear  energy  above  two  different  heights  as  a  function 
of  time  for  all  three  tubes.  This  normalized  energy  is  given  by 


ry=90Mm  rx=45Mm 

„  _  J y=yi  Jx=-45Mm  0zUXay 

^shear  rv=90 Mm  rx=45 Mm  ,  , 


(20) 


The  first  height,  iq  =  1.2  Mm,  is  the  lowest  height  at  which  /3  =  0.2  in  the  simulations  in  this  paper  with  initial 
non-zero  coronal  field,  and  is  an  estimate  of  the  height  at  which  we  need  to  transfer  a  significant  amount  of  the 
shear  magnetic  energy  relative  to  the  total  magnetic  energy  in  order  to  drive  a  CME.  The  second  height  is  at 
yi  =  5  Mm  above  the  surface  and  is  at  a  location  just  below  where  the  isolated  flux  rope  is  formed  from  the 
center  of  the  low  tension  tube. 


As  the  low  tension  case  develops  (t=1700  s  to  t=2200  s),  the  crescent  shaped  current  sheet  is  formed  and 
mass  and  shear  field  drain  to  the  ends,  below  1.2  Mm.  Thus  the  shear  energy  above  1.2  Mm  drops  as  shown  by 
the  green/dashed  curve  in  Figure  [5|b).  The  coalescence  of  plasmoids,  which  leads  to  a  concentration  of  shear 
flux,  then  creates  an  increase  in  Eshear  at  t=2300  s.  The  expulsion  of  these  plasmoids  along  the  current  sheet, 
down  below  1.2  Mm,  creates  the  sharp  decrease  seen  shortly  after  at  f=2400  s.  As  the  flux  rope  formed  from 
this  current  sheet  then  rises  into  the  corona,  carrying  shear  field  with  it,  the  shear  energy  above  5  Mm  rises,  as 
can  be  seen  in  the  green/dashed  plot  in  Figurc[9](c).  The  gradual  rise  in  Eshear  for  the  high  tension  cases  above 
1.2  Mm  after  t=2800  s  (red/dot-dashed  and  blue/solid  plots  in  Figure[9lb))  is  caused  by  reconnection  at  nearly 
vertical  current  sheets  which  form  at  the  edges  of  the  active  region.  This  reconnection  creates  plasmoids  which 
are  able  to  transfer  a  small  amount  of  shear  above  1.2  Mm  but  which  do  not  rise  to  the  5  Mm  level. 


The  magnetic  shear  energy  above  1.2  Mm  for  all  three  tubes  does  not  exceed  10%  of  the  total  magnetic 
energy  above  1.2  Mm.  For  the  emergence  of  a  2.5D  flux  tube  into  a  field-free  corona,  the  shear  field  is  generally 
confined  to  the  high  /3  lower  atmosphere.  The  shear  field  which  does  rise  in  the  coronal  flux  rope  is  not  strong 
enough  to  create  magnetic  pressure  which  will  counterbalance  the  magnetic  tension  and  gravitational  force  in  the 
rope  and  hence  drive  a  further  rise  or  an  eruption.  In  the  next  section  we  look  at  the  expansion  into  a  quadrupole 
field,  to  see  if  the  small  amount  of  shear  field  which  does  emerge  into  the  corona  can  drive  reconnection  at  the 
null  point  of  the  quadrupole,  thus  reconfiguring  the  field  to  allow  an  eruption  as  in  the  breakout  model. 


3.3.  Emergence  into  a  Quadrupole  Coronal  Field 

We  present  results  of  the  emergence  of  the  low  tension  flux  tube  into  a  pre-existing  coronal  quadrupole 
magnetic  field.  The  quadrupole  configuration  can  be  destabilized  by  the  presence  of  magnetic  shear,  as  in  the 
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breakout  model,  but  the  effects  of  the  lower  atmosphere  have  not  yet  been  investigated.  We  know  from  the 
previous  section  that  the  emergence  of  magnetic  flux  can  bring  shear  into  the  corona,  but  for  a  field-free  corona, 
the  coronal  flux  ropes  which  then  form  are  not  sheared  enough  to  erupt.  In  this  section  we  emerge  the  low 
tension  flux  tube  into  a  quadrupole  of  varying  strength,  and  examine  the  interaction  between  emerging  field  and 
coronal  field.  We  also  look  at  the  effect  of  the  quadrupole  strength  on  the  overall  emergence.  We  test  whether 
the  quadrupole  configuration  can  be  destabilized  by  flux  emergence-supplied  magnetic  shear  energy.  The  initial 
configuration  is  shown  in  Figure  [5] for  the  medium  strength  quadrupole  (B\x—Xa  /2,y=0  —  98  G). 

The  initial  rise  of  the  flux  tube  is  similar  to  the  field-free  corona  case  up  to  about  1000  s.  The  presence  of 
quadrupole  field  below  the  surface  has  a  negligible  effect  on  the  rise  of  the  tube,  and  its  initial  expansion  into  the 
atmosphere.  Figure  [10]  shows  the  evolution  of  the  low  tension  flux  tube  and  the  medium  strength  quadrupole 
at  6  different  times  (the  coupling  of  Bz  and  p  in  these  2.5D  simulations  allows  us  to  only  show  the  density  to 
demonstrate  the  structure  of  the  coronal  flux  ropes).  The  field  of  the  inner  dipole  of  the  quadrupole  is  aligned 
with  the  twist  field  of  the  upper  part  of  the  tube  which  expands  into  the  corona.  Thus  the  inner  flux  system 
of  the  quadrupole  field  is  pushed  upwards  towards  the  null  by  the  emergence,  as  shown  in  Figures  HOlal  and 
(b).  The  expanding  shell  of  the  tube  causes  a  crescent  shaped  current  sheet  to  be  formed  as  the  original  null  is 
stretched  over  this  shell.  Reconnection  sets  in  on  either  side  of  the  apex  of  the  current  sheet.  The  current  sheet 
continues  to  rise  as  the  expanding  shell  of  the  emerging  flux  tube  pushes  it  up  as  shown  in  Figures  [TOlVlWf'l. 
while  the  reconnection  creates  a  coronal  flux  rope  at  a  height  of  18  Mm. 

In  Figure  [10]  we  can  also  see  the  deformation  of  the  emerging  flux  tube’s  center,  and  the  creation  of  a  flux 
rope  by  reconnection  at  the  current  sheet  formed  at  about  5  Mm,  just  as  in  the  field-free  corona  simulations. 
Thus  there  are  two  coronal  flux  ropes  formed,  as  can  be  seen  in  Figure  ITOff).  The  upper  rope  is  created  by 
the  reconnection  at  the  separatrix  between  the  two  systems  of  the  quadrupole,  driven  by  the  expansion  of  the 
flux  tube  into  the  atmosphere.  The  lower  rope  is  formed  by  reconnection  at  the  current  sheet  formed  from  the 
center  of  the  emerging  flux  tube.  Both  of  these  flux  ropes  carry  shear  field  and  mass. 

Figure  |TT]  shows  the  vertical  forces  in  these  two  coronal  flux  ropes  at  later  times  of  t  =  2392  s  and 
t  =  2760  s.  Both  ropes  carry  significant  amounts  of  mass  and  in  both  ropes  the  magnetic  tension  acts  inwards, 
counteracting  any  outward  force  due  to  gas  or  magnetic  pressure.  Figure  [TTTal  shows  that  both  the  magnetic 
tension  and  magnetic  pressure  forces  are  directed  inwards  in  the  upper  flux  rope,  thus  a  strong  outward  pressure 
force  is  required  for  force  balance.  In  addition,  the  downward  gravity  force  of  the  trapped  mass  is  as  strong 
as  any  of  the  pressure  or  tension  forces.  These  other  forces  all  have  an  upward  directed  asymmetry  to  help 
balance  this  large  gravity  force,  but  nonetheless,  as  can  be  seen  in  Figure  [ITT  cl .  they  cannot  and  the  flux  rope 
starts  to  sink.  In  contrast,  the  gravity  force  is  relatively  small  in  the  lower  flux  rope,  as  shown  iiillllbl.  and  the 
magnetic  pressure  force  is  directed  outwards.  Yet  this  outward  magnetic  pressure  force  is  well  balanced  by  the 
inward  magnetic  tension.  Thus,  while  the  ropes  differ  in  the  way  the  magnetic  and  gas  pressures  act,  the  key 
fact  is  that  the  amount  of  shear  field  in  them  is  not  enough  to  create  an  unbalanced  outward  acting  magnetic 
pressure  force  to  cause  any  further  expansion. 

The  fraction  of  the  magnetic  energy  contained  in  the  shear  field  above  the  heights  of  1.2  Mm  and  5  Mm  for 
the  three  choices  of  quadrupole  and  the  field- free  corona  case  are  shown  in  Figure  IT21  In  general,  increasing  the 
strength  of  the  quadrupole  decreases  this  fraction  as  the  total  magnetic  energy  in  the  system  increases.  However 
this  also  increases  the  amount  of  reconnection  at  the  separatrix  of  the  quadrupole  as  the  configuration  contains 
more  flux  above  and  below  the  null. 

We  have  shown  that  the  emergence  of  a  flux  tube  into  a  quadrupole  field  does  cause  reconnection  at  the 
separatrix  between  the  two  flux  systems  in  the  quadrupole.  However,  the  nature  of  the  reconnection  means 
that  it  does  not  remove  significant  amounts  of  horizontal  field  from  above  and  below  the  original  null  point, 
as  in  the  single  X-point  reconnection  seen  in  the  breakout  model.  The  reconnection  at  the  null  is  driven  here 
by  expansion  of  the  outer  part  of  the  emerging  flux  tube,  and  the  current  sheet  formed  prefers  reconnection 
off-center,  creating  a  coronal  flux  rope  which  does  not  have  enough  magnetic  pressure  to  escape.  This  flux  rope 
eventually  falls  back  down,  and  joins  with  the  flux  rope  which  was  created  from  the  center  of  the  emerging  low 
tension  flux  tube.  The  magnetic  energy  in  the  shear  field  above  the  /3  =  0.2  level  is  too  small  (less  than  10%  of 
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the  total  magnetic  energy)  to  provide  the  shear  energy  required  by  the  breakout  model,  which  is  consistent  with 
the  fact  that  the  shear  field  in  the  flux  ropes  is  not  strong  enough  to  overcome  the  tension  and  gravitational 
forces. 

In  the  next  section  we  show  the  emergence  of  flux  tubes  into  a  dipole  coronal  field  which  is  orientated 
opposite  to  the  upper  part  of  the  emerging  flux  tube.  Reconnection  between  the  expanding  tube  and  the  pre¬ 
existing  dipole  coronal  field  may  remove  twist  field  in  the  coronal  flux  ropes  formed  and  allow  them  to  rise 
further. 


3.4.  Emergence  into  a  Dipole  Coronal  Field:  Effect  of  Tension  Profile 


In  this  section  we  present  results  from  simulations  of  the  emergence  of  the  three  different  flux  tubes  into  a 
particular  dipole  field  to  highlight  the  importance  of  the  choice  of  tube  profile.  We  choose  the  dipole  field  with 
B\x=o,y=o  =  52  G.  This  choice  of  dipole  best  highlights  the  differences  between  the  high  and  low  tension  cases. 
We  will  investigate  the  choice  of  dipole  strength  in 


Figure  fl3l  shows  the  emergence  of  the  high  and  low  tension  tubes  into  the  dipole  field.  The  behavior  of  the 
flux  tubes  prior  to  emergence  is  very  similar  to  the  field-free  corona  and  quadrupole  coronal  field  cases.  The 
outer  part  of  the  flux  tube  and  the  coronal  dipole  field  he  on  either  side  of  a  separatrix.  As  the  flux  tubes  rise 
into  the  dipole  field  this  separatrix  is  pushed  upwards  and  becomes  a  current  sheet,  as  shown  in  Figures Qj3a) 
and  (d).  In  both  cases  the  current  sheet  undergoes  reconnection,  consistent  with  the  tearing  mode  instability, 
creating  plasmoids  on  either  side  of  the  apex  which  are  expelled  along  the  current  sheet.  Eventually  a  coronal 
flux  rope  is  formed  by  coalescence  of  plasmoids  near  the  apex  of  the  current  sheet.  This  flux  rope  has  an  initial 
rise  speed  of  25  km/s.  In  the  low  tension  case,  the  flux  rope  is  formed  more  rapidly  and  at  a  higher  point  in  the 
corona.  Also,  a  flux  rope  is  formed  in  the  lower  corona  by  the  deformation  of  the  center  of  the  emerging  low 
tension  flux  tube,  leaving  two  flux  ropes  in  the  corona.  The  structure  of  these  two  flux  ropes  is  similar  to  that 
of  the  flux  ropes  formed  in  the  field-free  and  quadrupole  coronal  field  experiment,  with  magnetic  tension  and 
a  gravitational  force  acting  to  prevent  any  further  rise  or  expansion.  The  creation  of  these  coronal  flux  ropes 
appears  to  be  the  tearing  mode  instability,  as  was  the  case  in  simulations  of  flux  emergence  into  a  pre-existing 
field  performed  by  Archontis  et  al.  (2006t  20QZ)-  However  this  is  a  different  mechanism  from  that  which  creates 
the  flux  ropes  in  the  simulations  of  Manchester  et  al.  ( 2004 ) ,  Archontis  fe  Torok  ( 2008h ,  and  Fan  ( 2009 ) ,  which 
is  based  on  shear  flows,  converging  motions,  and  internal  reconnection  of  an  emerging  flux  tube. 

FigurefTHaf  shows  the  height  of  the  flux  tube  centers  (solid  lines)  and  the  height  of  the  separatrices  between 
flux  tube  and  dipole  (dashed  lines),  along  the  x  =  0  line  for  the  three  tubes.  Initially  the  paths  of  the  centers 
and  separatrices  are  the  same  for  all  three  tubes  as  they  expand  and  push  the  dipole  field  upwards.  After 
t  =  1500  s  the  low  tension  tube  undergoes  reconnection  at  its  center,  forming  the  lower  of  the  two  flux  ropes. 
Reconnection  at  the  separatrix  between  the  emerging  tube  and  the  dipole  field  in  all  three  simulations  creates 
the  coronal  rope  higher  in  the  corona  (solid  lines).  These  coronal  flux  ropes  eventually  fall  back  down  through 
the  transition  region  after  t  =  4000  s. 

The  amount  of  normalized  magnetic  shear  energy  supplied  to  the  two  heights  of  1.2  Mm  and  5  Mm  for 
all  three  tubes  emerging  into  the  dipole  can  be  seen  in  Figures  ITlTbl  and[FjJc).  The  normalized  shear  energy 
above  1.2  Mm  increases  after  2000  s  for  the  experiment  where  the  low  tension  tube  emerges  into  the  dipole, 
shown  as  the  green/dashed  curves  on  Figure  [mb').  This  is  caused  by  the  removal  of  some  horizontal  flux  by 
reconnection  in  the  current  sheets.  While  this  is  similar  to  the  reconnection  which  leads  to  an  eruption  in  the 
breakout  model,  this  reconnection  does  not  remove  enough  horizontal  field  to  raise  the  importance  of  the  shear 
field  in  the  coronal  flux  ropes  formed  by  this  reconnection.  As  a  result  the  shear  magnetic  energy  is  still  less 
than  10%  of  the  total  magnetic  energy  above  the  (3  =  0.2  corona.  3D  simulations  by  Archontis  fc  Torok  (20081) 
showed  that  reconnection  between  emerging  flux  tubes  and  horizontal  coronal  field  allows  coronal  flux  ropes 
which  are  formed  during  flux  emergence  to  escape.  This  is  not  the  case  in  these  2.5D  simulations  with  a  dipole 
coronal  field.  In  the  next  section  we  investigate  the  importance  of  the  dipole  strength  on  the  behavior  of  the 
coronal  flux  ropes. 
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3.5.  Emergence  into  a  Dipole  Coronal  Field:  Effect  of  dipole  field 

By  varying  the  dipole  strength  we  can  control  the  location  and  amount  of  reconnection  at  the  current 
sheet  formed  between  the  flux  tube  and  the  dipole,  and  through  this  affect  the  amount  of  shear  energy  (relative 
to  total  magnetic  energy)  which  is  transferred  to  the  low  j3  corona.  In  this  section  we  present  results  of  the 
emergence  of  the  low  tension  flux  tube  into  a  range  of  dipoles.  These  dipoles  have  horizontal  fluxes  which  are 
factors  of  the  horizontal  flux  that  is  initially  in  the  top  half  of  the  emerging  tube  of  1,  0.5  and  0.25.  We  also 
compare  these  three  cases  to  the  field-free  corona  case. 

The  strength  of  the  dipole  affects  the  expansion  of  the  outer  part  of  the  flux  tube,  as  the  stronger  the  dipole 
is  the  more  tension  it  has  to  resist  the  expansion  of  the  upper  part  of  the  flux  tube  into  the  corona.  The  stronger 
the  dipole,  the  quicker  (and  lower  down  in  the  corona)  the  two  systems  are  forced  together  by  the  emergence 
and  so  the  quicker  the  reconnection  occurs  between  them  and  the  quicker  the  coronal  flux  rope  is  formed.  This 
can  be  seen  in  Figure  [15]  which  shows  snapshots  at  t=2530  s  for  the  low  tension  flux  tube  emerging  into  the 
3  different  dipoles  (of  decreasing  strength  from  left  to  right)  and  a  field- free  corona  (far  right  panel).  At  this 
time,  the  stronger  dipole  case  has  already  formed  a  coronal  flux  rope,  while  the  weaker  dipole  case  has  only  just 
started  reconnection  between  the  emerging  flux  tube  and  the  dipole. 

FigurelTCl'a)  shows  the  heights  of  the  flux  tube  centers  and  the  flux  tube/dipole  separatrices  for  the  different 
choices  in  dipole  strength.  The  separatrices  all  reach  similar  heights,  but  do  so  at  different  times,  due  to  the 
differing  tension  in  the  dipoles  which  resists  the  pushing  up  of  the  dipole  by  the  expanding  flux  tube. 

Figures  Mb)  and  Me)  show  the  normalized  shear  energy  above  1.2  Mm  and  5  Mm.  In  general,  the 
magnetic  shear  energy  due  to  flux  emergence  is  confined  too  low  in  the  atmosphere  to  have  a  strong  enough 
effect  on  the  magnetic  field  in  the  corona.  The  shear  energy  is  still  always  below  10%  of  the  total  energy  above 
the  low  /?  corona.  Any  shear  field  that  is  transported  to  the  corona  by  coronal  flux  ropes  is  too  small  to  affect 
the  twist  field  at  these  heights.  The  coronal  flux  ropes  formed  are  not  sheared  enough  and  contain  significant 
amounts  of  dense  material.  In  order  to  drive  a  CME  by  flux  emergence  in  these  simulations,  we  need  a  method 
for  transferring  more  shear  field  to  coronal  heights  while  transferring  less  mass. 


4.  CONCLUSIONS 


We  have  investigated  the  initiation  of  CMEs  by  dynamical  flux  emergence  using  2.5D  cartesian  simulations 
of  the  emergence  of  a  range  of  sub-surface  flux  tubes  into  a  range  of  coronal  magnetic  field  configurations.  Most 
current  CME  models  do  not  include  the  lower  dense  atmosphere  in  their  models,  but  the  boundary  conditions 
for  the  shear  energy  and  dynamical  flows  they  use  are  most  likely  results  of  flux  emergence  through  this  region. 
In  a  first  attempt  to  couple  the  modeling  of  CME  initiation  and  flux  emergence,  we  have  taken  coronal  field 
configurations  which  are  based  on  the  magnetic  breakout  model.  However,  our  results  are  generally  applicable 
to  any  CME  model  which  assumes  the  appearance  of  highly  sheared  flux  at  the  low  /3  corona. 


The  evolution  of  flux  tubes  into  a  field-free  corona  can  be  separated  into  low  and  high  tension  cases,  as  in 
Murray  fc  Hood  (2008il.  based  on  the  relative  amount  of  tension  near  the  center  of  the  flux  tube.  During  the 


emergence  of  a  low  tension  tube,  an  isolated  flux  rope  was  created  by  the  deformation  of  the  flux  tube  center  into 
a  current  sheet  and  the  reconnection  of  fieldlines  on  either  side  of  the  apex  of  this  current  sheet.  Since  twisted 
flux  ropes  are  essential  for  several  CME  models  dTitov  fe  Demoulinlll999t  iRoussev  et  alJl2003t  iKliem  fe  Torok 


20061:  Fan  &  Gibson  20071)  further  investigations  are  warranted,  especially  in  3D,  to  take  into  account  the  effect 


of  curvature  along  the  tube’s  axis.  The  comparison  of  3D  simulations  with  observations  of  neutral  lines  during 
flux  emergence  will  be  able  to  verify  if  this  behavior  is  occurring  on  the  Sun. 


In  an  attempt  to  test  whether  the  shear  field  brought  into  the  corona  by  flux  emergence  could  destabilize 
the  quadrupole  configuration  of  the  magnetic  breakout  model,  we  added  a  pre-existing  quadrupole  coronal  field 
above  the  emerging  flux  tubes.  The  expanding  outer  shell  of  the  flux  tube  stretched  the  original  null  point  of  the 
quadrupole  into  a  current  sheet,  and  reconnection  in  this  current  sheet  created  another  coronal  flux  rope,  higher 
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in  the  corona.  Both  types  of  flux  ropes  were  created  in  similar  fashions:  by  the  formation  of  a  crescent  shaped 
current  sheet,  with  preferential  reconnection  occurring  off-center,  forming  a  coronal  flux  rope  at  the  apex  of  the 
crescent.  However,  one  type  of  flux  rope  was  created  by  the  deformation  of  the  center  of  an  emerging  flux  tube, 
and  one  was  created  along  the  separatrix  between  the  two  flux  systems  of  the  quadrupole. 

The  emergence  of  sub-surface  flux  tubes  into  a  simpler  dipole  coronal  field  configuration  was  simulated, 
in  an  attempt  to  remove  the  twist  field  of  the  flux  tube  which  first  emerges  into  the  corona.  Reconnection 
at  the  separatrix  between  emerging  field  and  pre-existing  dipole  field  yielded  an  isolated  flux  rope,  as  did  the 
deformation  of  the  center  of  the  emerging  low  tension  flux  tube 

Although  the  coronal  flux  ropes  that  are  formed  in  these  simulations  of  flux  emergence  transport  shear 
field  into  the  low  (3  corona,  the  shear  field  is  weak  compared  to  the  twist  field  and  they  are  laden  with  photo- 
spheric/chromospheric  material.  Any  outward  directed  forces  from  the  shear  field  are  insufficient  to  overcome 
the  inward  directed  tension  force  and  the  downward  directed  gravity  force.  These  coronal  flux  ropes  are  therefore 
not  candidates  for  erupting  structures,  and  eventually  disappear  or  fall  back  into  the  lower  atmosphere. 

We  have  covered  a  range  of  flux  tube  profiles  and  coronal  field  structures  and  strengths,  and  in  all  cases 
found  that  the  emergence  of  sub-surface  flux  tubes  in  these  2.5D  simulations  is  an  inefficient  method  for  the 
transfer  of  magnetic  shear  energy  into  the  low  (3  corona.  The  amount  of  magnetic  shear  energy  supplied  to  the 
low  (3  corona  was  calculated,  at  heights  of  1.2  Mm  and  5  Mm.  The  lower  of  these  two  heights  corresponds  to  the 
lowest  point  at  which  (3  =  0.2  in  the  simulations  in  this  paper,  and  thus  the  lowest  point  at  which  the  conditions 
are  the  same  as  those  used  in  CME  models.  In  all  cases,  the  amount  of  magnetic  shear  energy  supplied  to  the 
/ 3  <  0.2  corona  was  less  than  10%  of  the  total  magnetic  energy  and  appears  insufficient  to  directly  drive  a  CME. 
We  conclude  that  the  simple  observation  of  flux  emergence  and  shear  at  the  photosphere  is  not  sufficient  evidence 
that  a  CME  can  be  driven.  Further  investigations  must  be  carried  out  to  help  understand  the  conditions  under 
which  emerging  shear  can  be  transported  up  to  the  low  corona,  where  it  is  needed. 

For  3D  flux  tubes  with  significant  curvature  along  the  tube  axis,  flows  along  fieldlines  may  be  an  important 
phenomena  in  flux  emergence  simulations.  These  simulations  should  be  repeated  in  3D  to  see  if  this  has  an 
affect  on  the  formation  of  coronal  flux  ropes  and  their  eventual  evolution  in  the  solar  atmosphere. 


We  have  tested  the  robustness  of  these  results  by  varying  the  resolution  and  introducing  asymmetry  into 
the  flux  emergence  and  coronal  fields,  and  find  similar  results  regardless  of  the  position  of  the  flux  tube  relative 
to  the  center  of  the  quadrupole/dipole. 


In  these  simulations,  and  almost  all  previous  flux  emergence  simulations,  the  high  (3  lower  atmosphere 
is  assumed  to  be  fully  ionized,  which  is  not  the  case  for  the  solar  chromosphere.  The  presence  of  neutrals 
in  the  nearly  neutral  region  near  the  temperature  minimum  may  have  an  important  effect  on  these  results. 
Leake  fe  Arber  ( 2006h  and  Arber  et  al.  (2007)  showed  that  this  nearly  neutral  region  can  dramatically  affect 


flux  emergence.  In  particular,  cross-field  currents  in  emerging  structures  are  destroyed  by  ion-neutral  collisions, 
and  the  magnetic  field  lines  are  able  to  slip  through  the  neutral  gas  during  emergence.  Thus  the  flux  tubes  lift 
up  less  mass  from  the  chromosphere  when  the  effects  of  partial  ionization  are  taken  into  account.  All  this  occurs 
lower  down  in  the  solar  atmosphere  around  the  (3=1  region  and  will  directly  affect  the  transfer  of  both  mass 
and  shear  energy  to  the  low  (3  corona  above.  These  effects,  along  with  3D  effects,  may  provide  a  mechanism  for 
the  transfer  of  magnetic  shear  energy  into  the  corona  without  the  transfer  of  mass,  which  may  be  key  to  driving 
a  CME  with  dynamical  flux  emergence.  This  topic  will  be  investigated  in  future  publications. 
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Fig.  1. —  Meridional  projections  of  magnetic  field  lines  throughout  the  magnetic  breakout  eruption  process  at 
six  different  times.  Reprinted  with  permission  from  Lynch  ed  ah  (120041 .  Spatial  units  are  solar  radii. 


19  - 


- -P9 

5  1-5 

(/) 

- “P9 

_ +B-VB/47T 

C  1.0 

>N 

7/  \  _ +B-VB/47T 

_ -V(B2)/8tt 

"O 

1  \  - -V(B2)/8^t 

o 

- 

'  0.5 

-VP 

O 

\  -VP 

_ 

"* — ' 

■  -  - - ...  

<D  0.0 
o 

o 

^  •  ■ — .... _ 

time  =  0.00000  s 

time  =  30000.0  s 

1.05  1.10  1.15  1 .20  ^  1.05  1.10  1.15  1.20 


Solor  Rodii  Solor  Rodii 


0  20  40  60  80  100 

Time  (1000  s) 


Fig.  2.-  Panels  (a)  and  (b):  Vertical  forces  along  the  latitude=0°  line  at  the  center  of  the  inner  arcade  of  the 
magnetic  quadrupole  in  the  breakout  model  at  t=0  s,  panel  (a),  and  t=30000  s,  panel  (b),  after  shearing  of 
the  arcade  by  imposed  boundary  flows.  Panel  (c):  The  normalized  shear  magnetic  energy  (solid  line)  and  2D 
integral  of  kinetic  energy  density  (dashed  line)  in  the  system  during  this  shearing  period. 
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Fig.  3.-  Panel  (a):  the  initial  height  profile  of  the  model  stratified  atmosphere  showing  the  temperature  (solid 
line),  gas  pressure  (dashed  line),  and  magnetic  pressure  (dot-dashed  line),  up  to  20  Mm.  The  magnetic  pressure 
profile  is  for  a  sub-surface  magnetic  flux  tube  and  an  overlying,  anti- parallel  dipole.  Panel  (b):  The  initial 
plasma  f3  for  six  different  simulations:  Three  simulations  with  background  dipole  field  and  three  simulations 
with  background  quadrupole  field.  The  straight  black  dotted  line  shows  a  constant  value  of  /3  =  0.2. 
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Fig.  4.—  Panel  (a):  Initial  twist  function,  0,  for  the  three  flux  tubes.  Tube  1:  Constant  twist,  high  tension, 
blue/solid  lines.  Tube  2:  Increasing  twist,  low  tension,  green/dashed  lines.  Tube  3:  Decreasing  twist,  high 
tension,  red/do t-dashed  lines.  Panel  (b):  The  axial  field,  BZl  and  twist  held,  Bg ,  for  the  same  three  tubes. 
Panel  (c):  The  gradient  in  magnetic  pressure,  —V(B2)/8tt  (dashed  lines),  and  the  magnetic  tension,  B.\7B/1-k 
(solid  lines),  both  normalized  by  Bq/cl ,  for  the  same  three  tubes. 
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Fig.  5. —  Log  of  density  at  t  =  0  and  initial  magnetic  field  lines  in  the  x  —  y  plane  (constant  contour  intervals 
of  the  z  component  of  the  magnetic  vector  potential,  Az).  The  field  lines  show  the  outer  part  of  the  sub-surface 
flux  tube  (located  at  x=0,  y=-2  Mm),  and  an  overlying  quadrupole  field.  The  inner  arcade  of  the  quadrupole 
is  parallel  to  the  upper  half  of  the  flux  tube,  as  can  be  seen  by  the  arrows,  which  show  the  direction  of  Bx  in 
the  plane.  The  null  is  located  at  about  3.2  Mm  above  the  surface  (y=0). 


-  23  - 


Fig.  6.  Panels  (a)  to  (d):  The  magnetic  shear  field,  Bz  (G),  and  magnetic  held  lines  (constant  contour  intervals 
of  Az)  at  four  different  times  for  the  high  tension  flux  tube  (Tube  1)  emerging  into  a  field-free  corona.  Panels 
(e)  to  (h):  The  same  but  for  the  low  tension  tube  (Tube  2). 
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Fig.  7.  —  Panels  (a)  to  (c):  The  magnetic  shear  field,  Bz  (G),  and  magnetic  field  lines,  at  three  different  times 
for  the  low  tension  flux  tube  emerging  into  a  field-free  corona.  Panels  (d)  to  (f):  The  same  but  for  the  log  of 
the  density.  The  two  white  tick  marks  on  the  left  are  located  at  y=1.2  Mm  and  5  Mm. 
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Fig.  8.  Panels  (a)  and  (b):  Vertical  forces  along  the  central,  x  =  0,  line  inside  the  initial  buoyant  sub-surface 
flux  tubes  at  t  =  0  s  for  the  high  tension  and  low  tension  cases.  Panels  (c)  and  (d):  Vertical  forces  along  the 
central,  x  =  0,  line  above  the  surface  at  t  =  2760  s  for  the  high  tension  and  low  tension  cases.  Panel  (d)  shows 
the  coronal  flux  rope  formed  from  the  emerging  low  tension  tube. 
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Fig.  9. —  Panel  (a):  The  height  of  the  center  of  the  three  flux  tubes  emerging  into  a  field-free  corona.  Panels  (b) 
and  (c):  The  normalized  magnetic  shear  energy,  -Eghegj  (see  Eauationl2Ul) .  above  1.2  Mm  and  5  Mm  respectively, 
for  the  same  three  tubes. 
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Fig.  10. —  Log  of  density  and  fieldlines  for  the  low  tension  tube  at  six  different  times  as  it  emerges  into  a 
medium  strength  (B \x—Xa/2,y=o  =  98  G)  quadrupole  field.  The  arrows  show  the  direction  of  Bx  in  the  x  —  y 
plane. 
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Fig.  11. —  Vertical  forces  along  the  central,  x  =  0,  line  for  the  two  flux  ropes  formed  when  the  low  tension 
flux  tube  emerges  into  the  medium  strength  {B\x—Xa /2,y=o  =  98  G)  quadrupole.  Panel  (a)  shows  the  forces  in 
the  upper  flux  rope  formed  by  reconnection  at  the  current  sheet  formed  from  the  initial  quadrupole  null  point. 
Panel  (b)  shows  the  forces  for  the  flux  rope  created  by  the  deformation  of  the  center  of  the  emerging  flux  tube. 
Panels  (c)  and  (d)  show  the  total  vertical  forces  in  these  two  flux  ropes,  respectvely. 
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Fig.  12. —  Normalized  magnetic  shear  energy,  Egjiear  (see  Equation [20li .  above  two  different  heights  for  the  low 
tension  flux  tube  emerging  into  a  quadrupolar  coronal  field.  Panel  (a)  is  for  1.2  Mm.  Panel  (b)  is  for  5  Mm. 
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Fig.  13.'  Log  of  density  and  fieldlines  at  three  different  times  for  the  emergence  of  a  high  tension  flux  tube, 
panels  (a)  to  (c),  and  a  low  tension  flux  tube,  panels  (d)  to  (f),  as  they  emerge  into  the  medium  strength, 
B\x-q =  52  G,  dipole  coronal  field.  The  arrows  show  the  direction  of  the  magnetic  field. 
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Fig.  14.—  Panel  (a):  Height  of  the  flux  tube  centers  (dotted  lines)  and  tube/dipole  separatrices  (solid  lines) 
for  the  three  flux  tubes  emerging  into  a  dipole  of  medium  strength  (B \x=0  y=0  =  52  G).  Panels(b)  and  (c):  The 
normalized  magnetic  shear  energy,  Egjlear  (see  Equation  l20l).  above  heights  1.2  Mm  and  5  Mm  for  the  same 
three  experiments. 
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Fig.  15. —  Shear  field,  Bz{ G),  and  fieldlines  at  t=2530  s  for  the  emergence  of  the  low  tension  flux  tube  into 
4  coronal  fields  of  decreasing  strength.  Panel  a):  B |x=o,y=o  =  104  G.  Panel  b):  B\x=o,y=o  =  52  G.  Panel  c): 
B\x=0,y=o  =  26  G.  Panel  d):  B \x=0  y=0  =  0  G  (no  dipole).  Panels  (e)  to  (h)  show  the  same  but  for  the  log  of 
density.  Arrows  show  direction  of  magnetic  field. 
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Fig.  16.-  Panel  (a):  Height  of  the  flux  tube  centers,  (dashed  lines)  and  tube/dipole  separatrices  (solid  lines) 
for  the  low  tension  flux  tube  emerging  into  dipoles  of  varying  strength.  Panels  (b)  and  (c):  The  normalized 
magnetic  shear  energy,  Es jlear  (see  Equation  EUl).  above  heights  1.2  Mm  and  5  Mm  for  the  same  experiments. 
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